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Abstract. Recent far-infrared and submillimetre waveband observations revealed large number of Ultraluminous 
Infrared Galaxies (ULIGs) with infrared luminosities > 10 12 Lq. These sources are proposed to lie at redshifts 
above one, and in normally interacting systems with very dusty environments. We discussed in a previous paper 
that a population with a fast evolving infrared burst phase triggered by gas-rich mergers at z ~ 1 interpreted 
successfully the steep slope of faint IRAS 60 pm source counts within the flux range of 100 mJy ~ 1 Jy, still 
leaving the infrared background level at this wavelength compatible with the upper limit from recent high energy 
TeV 7 ray detection of Mrk501. To extend the model to mid and far infrared wavelengths, we adopt a reasonable 
template spectral energy distribution as typical for nearby infrared bright starburst galaxies ( Li r <= 10 12 Lq), 
such as Arp220. We construct the SED for the dusty starburst mergers at z ~ 1 by a simple dust extinction law 
and a thermal continuum assumption for the far-infrared emission. Since the radiation process at mid-infrared 
for these starburst merging systems is still uncertain, we assume it is similar to the MIR continuum of Arp220, 
but modify it by the observed flux correlation of ULIGs from IRAS and ISOCAM deep surveys. We show in this 
paper that the strong evolution of the European Large Area ISO Survey (ELAIS) at 90 pm, ISO 170 pm and 
the Submillimeter deep survey at 850 pm could be sufficiently accounted for by such an evolutionary scenario, 
especially the hump of the ISOCAM 15 pm source count around OAmJy. From current best fit results, we find 
that the dust temperature of those extremely bright starburst merging system at z ~ 1 would be higher than that 
of Arp220 for a reconciliation of the multi-wavelength infrared deep surveys. We thus propose that the infrared 
burst phase of dusty starburst galaxies or AGNs from gas-rich mergers at z ~ 1 could contribute significantly to 
the strong evolution of the IRAS 60 pm, the ISO 15 pm, 90 pm, 170 pm, as well as the SCUBA 850 pm number 
counts, while being compatible with the current observational limits of cosmic infrared background and the 
redshift distributions. The major difference of our current model prediction is that we see a fast convergence of 
the differential number counts at 60 pm below 50 mJy, which is about a factor of two brighter than other model 
predictions. Future infrared satellites like Astro-F or SIRTF would give strong constraints to the models. 
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1. Introduction 

There has been much progress in the study of extragalactic 
evolution since far-infrared and submillimeter deep sur¬ 
veys detected a significant population of Ultraluminous 
Infrared Galaxies (ULIGs: Li r > 10 12 L@) at high red- 
shift (z ~ 1 — 4) (Hughes et al. 1998, 2000, Blain et al. 
1999, Eales et al. 1999, Holland et al. 1999, Puget et al. 
1999, Sanders 1999, Dole et al. 2000). The major inter¬ 
ests of present research are the nature and evolution of 
these sources. Due to the lack of high resolution morpho¬ 
logical studies, the origin of these faint SCUBA sources 
are still not clearly understood. Local ULIGs most likely 
arise from major galaxy mergers with high dust extinc- 
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tion for the central starburst or AGN activities (Sanders 
& Mirabel 1996, Murphy et al. 1996, Surace et al. 1998, 
2000, Scoville et al. 2000). Nevertheless, more than 50% 
of the high redshift ultraluminous infrared objects and 
the optical counterparts of ISOCAM HDF-N galaxies are 
suggested to show indication of galactic interactions and 
merger signatures (Mann et al. 1997, Smail et al. 1999, 
Ivison et al. 1998, 2000, Sanders 1999). Meanwhile, esti¬ 
mation of their spectral energy distributions suggests that 
these galaxies are probably the high redshift counterpart 
of the local ULIGs discovered by the IRAS deep survey 
(Barger et al. 1998,1999, Smail et al. 1998, Trentham et al. 
1999, Frayer et al. 2000). Although the exact mechanism 
for formation of these interesting sources are not clear, 
there are certain reasons to say that majority of them may 
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be major mergers of gas-rich disks accompanied by dust- 
shrouded nuclear starbursts or powerful Active Galactic 
Nuclei (AGN). SMM J14011+0252, and ERO J164502- 
4626.4 (HR10) are two such candidates with their cen¬ 
tral activities heavily hidden by dust extinction, and are 
suggested to be consistent with the evolutionary track of 
mergers-starbursts/AGN, probably elliptical galaxies in 
formation (Graham & Dey 1996, Cimatti et al. 1998, 1999, 
Frayer et al. 1998, 1999, Dey et al. 1999, Papadopoulous 
et al. 2001 ). 

On the other hand, the source counts from present in¬ 
frared and submilimetre surveys, such as IRAS, ISO and 
SCUBA all significantly exceed the non-evolving predic¬ 
tions. The extremely strong evolution is seen from the 
differential counts of the ISOCAM at 15 /rm, with the re¬ 
markable upturn at £15 < 3 mJy and a fast convergence 
when £15 ~ 0.3 mJy. This striking feature is based on the 
data from several independent sky surveys (Elbaz et al. 
1999, Chary & Elbaz 2001, Mazzei et al. 2001, Serjeant et 
al. 2001). Although there are many other possible evolu¬ 
tionary scenarios which could explain the present obser¬ 
vations, the reason that we are encouraged to explore here 
a merger-driven galaxy evolution picture with binary ag¬ 
gregation dynamics is simply because the IRAS database 
and recent ISO and sub-mm deep surveys indicate that 
most of the luminous infrared sources are actually inter¬ 
acting/mergering systems. Also, the local IR luminosity 
function shows an excess over the Press-Schechter for¬ 
mula (Press & Schechter 1974, Lonsdale 1995, Pearson 
& Rowan-Robinson 1996, Guiderdoni et al. 1998, Roche 
et al. 1998, Rowan-Robinson et al. 1998, Dey et al. 1999, 
Sanders 1999, Dole et al. 2000, Efstathiou et al. 2000, Silk 
& Devriendt 2000, Serjeant et al. 2001, Takeuchi et al. 
2001 ). 

Considering mergers as a possible formation mecha¬ 
nism of Ultraluminous Infrared Galaxies both at high and 
low redshift, as well as their significant infrared emissions, 
Wang (1999) and Wang & Biermann (2000) discussed 
the effects of galaxy mergers on the strong evolution of 
the IRAS 60 ym deep survey within a binary aggregation 
galactic evolutionary scheme. In this model, the bright tail 
of the infrared luminosity function is simulated in a con¬ 
sistent way for both the density and luminosity evolution 
due to the decrease of the merger fration with cosmic time 
and a merger-triggered infrared burst phase. They found 
that a luminosity-dependent infrared burst phase is cru¬ 
cial for the interpretation of the steep slope within a flux 
range of 10 mJy ~ 1 Jy by the IRAS 60 ym deep sur¬ 
vey. This means dusty starburst galaxies or AGNs from 
gas rich mergers at high redshift may experience an in¬ 
frared burst phase around a transition redshift z ~ 1 , and 
fade quickly within the merger time scale of that epoch. 
The more massive merger systems could have such in¬ 
frared emission enhanced to a higher level and decrease 
even faster. This kind of speculation is based on the ob¬ 
servation that ULIGs are usually more than a factor of 20 
brighter than normal starburst galaxies. Although the de¬ 
tailed mechanism for such enormous infrared emission is 


still unclear, it is believed to be related to a special stage of 
the merger process when the dust mass and temperature 
are both dramatically increased (Kleinmann & Keel, 1987, 
Taniguchi & Ohyama 1998). Recent numerical simulation 
on the evolution of dusty starburst galaxies by Bekki & 
Shioya (2001) shows that there is a very strong photomet¬ 
ric evolution during the merger process of two gas-rich 
disks, and a dramatic change of the spectral energy dis¬ 
tribution (SED) around a cosmic time scale T ~ 1.3 Gyr, 
when the two disks of the merger become very close and 
suffer from violent relaxation and the star formation be¬ 
comes maximal (~ 378 Mgj/r^ 1 ). The infrared flux in this 
case could increase by one magnitude, especially for the 
far infrared wavelength range (60 ym ~ 90 ym) in the 
emitting frame. 

The redshift distribution of the contributing sources 
for the steep slope at faint IRAS 60 ym counts in the 
model of Wang & Biermann (2000) shows that the in¬ 
frared burst phase around z ~ 1 could have comparable 
significance to the local IR sources. The question is then 
whether such an infrared burst phase, or such a popu¬ 
lation of ULIGs, could also sufficiently account for the 
strong evolution seen in other infrared wavelengths, espe¬ 
cially at the ISOCAM 15 ym, ISOPHOT 90 ym, 170 ym 
and SCUBA 850 ym. We thus try to make a reconcilia- 
tory evolution model which could fit at least the present 
statistics of the multi-wavelength deep surveys. 

In this paper, we will first review the binary aggrega¬ 
tion galaxy evolution model by Wang (1999) and by Wang 
& Biermann (2000) in section 2, where starburst/AGN ac¬ 
tivities may be triggered during the merger process, as well 
as an infrared burst phase from gas-rich mergers around a 
redshift of one. We will discuss the SED template we adopt 
in our calculation for the nearby starburst galaxies and a 
possible strong evolution of the spectral energy distribu¬ 
tion of the dusty starburst merging systems at z ~ 1. We 
thus could further investigate whether the infrared burst 
phase from gas-rich mergers around redshift z ~ 1 is suffi¬ 
cient to account for the strong evolution also detected by 
ISO and submillimetre deep surveys. One set of cosmo¬ 
logical parameters, namely Hq = 50 km/s/Mpc, O = 0.3 
and A = 0.7 is adopted in the calculation. 

2. Model 

We adopt in this study the binary aggregation dynam¬ 
ics based on the Smoluchowski equation (1916) within a 
merger-driven galaxy evolutionary scheme, which simu¬ 
lates the evolution of a luminosity function due to galaxy 
mergers and predicts the redshift-dependent luminosity 
function of a population of galaxies evolving forwards in 
time from the formation epoch to match the observed 
local luminosity function, number counts and space dis¬ 
tribution etc. Studies of Cavaliere & Menci (1993, 1997) 
show that this method could include more dynamics to de¬ 
scribe a further step in galaxy-galaxy interactions within 
the scheme of direct hierarchical clustering (DHCs), and 
probably could help to alleviate some intrinsic problems 
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in the DHC scenario, such as the overproduction of small 
objects as well as the difficulty of a reconciliation be¬ 
tween the excess of faint blue counts and the flat local 
luminosity function. The numerical technique to solve the 
Smoluchowsky equation is a Monte-Carlo approach for 
the inverse-cascade merger tree. Readers are referred to 
Cavaliere & Menci (1993, 1997), Wang (1999) and Wang 
& Biermann (2000) for the details of the dynamics and 
the techniques. 

Considering different evolutionary characteristics of 
different morphologies, we adopt in our study a multi- 
component model that contains starburst galaxies, dust 
shrouded AGNs and spiral galaxies as three major classes 
of infrared emitting sources. The local luminosity func¬ 
tions of the spiral and starburst galaxies at 60 ym from 
Saunders (1990), and that of Seyferts from Ruch et al. 
(1993) are used to normalize the Monte-Carlo simulation. 
We adopt the mass-light relation of blue starburst galax¬ 
ies, which is given by Cavaliere & Menci (1997) in a study 
of the excess of faint blue galaxies in optical surveys. The 
abundances of dust-shrouded AGNs are set to be 50% and 
80% at local and high redshift based on the statistics from 
Hubble Space Telescope imaging survey of nearby AGNs 
and the cosmic X-ray background (Malkan et al. 1998, 
Gilli et al. 1999). 

The modelling of a luminosity-dependent ultralumi- 
nous infrared burst phase from gas-rich mergers is de¬ 
scribed in detail by Wang & Biermann (2000). Here we 
give a brief review of the basic dynamics and the tem¬ 
plate SED we adopt in our model for the infrared lumi¬ 
nous sources. We introduce also in this section the con¬ 
struction of a SED for the dusty starburst mergers at 
z ~ 1, normally with the luminosity Li r > 10 12 L@, in 
order to further investigate such an evolutionary scenario 
at mid and far-infrared wavelengths from the ISO deep 
survey. Considering that star formation is triggered by 
mergers and proportional to M gas /T (r is the dynam¬ 
ical interaction time scale), Cavaliere & Menci derived 
a mass-light ratio for dwarf galaxies, L/L* = (M/M*) v 
(where r/ = 4/3, if the cross section is purely geometri¬ 
cal). L* oc f(z, Ao,flo) could be used to describe a red- 
shift dimming, or a luminosity evolution. A power law 
prescription of f(z) oc (1 + z) 13 is adopted in the model. 
Simplifying the color and K- corrections, they roughly get 
Lb °c -py • We assume in our model a 

luminosity ratio oc M^, which is consistent with cur¬ 
rent understanding of the nature of ULIGs, where people 
normally believe that the extremely infrared bright phase 
is due to the starburst merger events with the far-infrared 
luminosity Li r enhanced both by the accumulation of the 
dust mass and the increase of the dust temperature. This 
burst phase could enhance the infrared luminosity by a 
factor of about 20 over that of normal starburst galax¬ 
ies (Kleinmann & Keel, 1987, Taniguchi & Ohyama 1998, 
Bekki & Shioya 2001). We give L eo oc f{z)M TlJril . 
y = 1 ~ 1.2 is adopted in the calculation, which not only 
reasonably represents an infrared enhancement of about a 


factor of 20 for a typical ULIG with a mass of 10 12 M 0 (the 
mass increases about one magnitude over that of normal 
starburst galaxies), but successfully interprets the steep 
slope of IRAS number counts. The scaling factor of the 
mass-light ratio here is normalized by the local luminosity 
function of the IRAS deep survey. We found from our best 
fit results that a population of infrared starburst sources, 
especially with spheroidal morphology, would have experi¬ 
enced very strong evolution in the past. A rate of (3 = 3.7 
in the luminosity evolution f(z) oc (1 + z)@ since a tran¬ 
sition redshift z ~ 1 indicates a very strong evolution for 
such a population of starburst galaxies. This is at least 
comparable to, if not stronger than, QSOs (Roche et al. 
1998, Lonsdale 1995, Pearson & Rowan-Robinson 1996, 
Rowan-Robinson et al. 1998, Sanders 1999, Dole et al. 
2000, Franceschini et al. 1988, 2001). A differential dim¬ 
ming is simulated by Li r (z — Sz) oc Li r (z ) 1_ ‘> below a 
transition redshift z ~ 1, in order to match the observed 
local luminosity function by the IRAS deep survey. The 
simulation gives a best value £ ~ 0.4. This power law sup¬ 
pression also includes another physical reality, that the in¬ 
frared luminous galaxies at the bright tail of the luminos¬ 
ity function become gas poor faster than the less luminous 
ones. Besides the merger rate decrease with cosmic time, 
this physical effect is very important for a good fit of the 
steep slope at IRAS 60 /im number counts within a flux 
range of 100 mJy ~ 1 Jy. 


We reviewed above the dynamics and some important 
physical parameters in the current study, which are the 
same as those used in the previous paper of fitting IRAS 
60 ym number counts (Wang 1999, Wang & Biermann 
2000). In the following, we will start to construct the 
spectral energy distribution for dusty starburst mergers 
around redshift z ~ 1, in order to extrapolate the calcu¬ 
lation to mid- and far-infrared wavelengths. Because the 
unclear nature of the Ultraluminous Infrared Galaxies at 
high redshift, we do not have a good understanding of 
the dust environment and properties in these sources. An 
optically thin, single-temperature dust model is adopted 
as a first order approximation for a modified blackbody 
continuum of temperature T at far-infrared wavelength 
in this calculation. The formula is simply given by S\ = 
B\(T)t\ oc B\K\. t\ = K\pdl is the dust opacity, 
and K\ the dust absorption coefficient ( K\ oc A _/ii of 
(3 ~ 1 — 2). In this case, the flux received at wavelength 

A is S\ = V ±‘ d B 2 ■ where T is the scaling factor for a 

conservation of the dust absorbed energy and re-emitting 

energy, A e = A/(1 + z) is the wavelength in the emitting 

frame. A flux ratio in the observer’s frame could be derived 

by ~ (^ ) d + 5 e ^(^-^)^ j where 

J S\ 2 A 2 /b(A 2 «,T) 

h is Plank’s constant, k is Boltzmann’s constant and C 
is the speed of light. With a reasonable assumption for 
the dust emissivity power (3 and the dust temperature T, 
we can easily extend our calculation to far-infrared wave¬ 
lengths. The mid-infrared emission is more complicated 
than that of the far-infrared which could be well described 
by a single temperature blackbody spectrum by cold, large 
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grain dust. The MIR emission properties are usually dom¬ 
inated by the radiation field of heated small grains and 
PAHs. These dust grains are normally heated stochasti¬ 
cally, and are not in thermal equilibrium with the ambi¬ 
ent radiation field. Thus the MIR continuum is most like a 
power law spectrum. In this calculation, we will not go to 
a detailed modelling of the MIR emission feature. Instead, 
we only modify the template SED of the starburst galaxy 
Arp220 by Silva et al. (1998), with the observational corre¬ 
lations of *Si 5 /*5*60 by IRAS and ISOCAM deep surveys for 
the ultraluminous case (> 10 11 Lq) to represent the 
dusty starburst merging system around redshift one. The 
color ratio of S 15 / Sqq for the ULIGs are a factor of about 
5 lower than the mean value of the whole sample, which 
may imply a very complicated process to heat small grains 
during the merger process (Aussel et al. 2000, Dunne 
et al. 2000, Saunders et al. 2000, Chary & Elbaz 2001). 
Although there are many indications that the MIR con¬ 
tinuum is correlated with the temperature of large grain 
dust (Dale et al. 2001), there is still no exact modelling 
for such a process. Our goal in this paper is to construct 
a simple SED for those luminous starburst mergers based 
on various obsevational correlations, which not only repre¬ 
sents the observed trend for individual samples of a certain 
luminosity bin {L ir > 10 12 Lq), but matches the statisti¬ 
cal results from the multiwavelength deep surveys. The 
number of sources dN in a comoving volume dV within 
the flux range S\ to S\ + dS\, measured at wavelength A, 
is defined by: dN = p(L x , z) dV ^ dS x , 

where K{L x ,z) = L is the K-correction, and dj, is 
the luminosity distance in a A dominated universe. 

3. Results and discussions 

The exact broad-band spectra of faint IR sources is still 
not well defined. Considering deep surveys at various 
IR/submm wavelengths would help to simultaneously con¬ 
strain the evolution properties and the typical spectral 
energy distribution of such sources. We show in this sec¬ 
tion the comparison of our model prediction with the 
ISOCAM 15 pm survey data, IRAS 60 pm, ISOPHOT 
90 pm and 170 pm (FIRBACK Survey), as well as the 
SCUBA 850 pm data. We also calculated the redshift dis¬ 
tribution of these sources within a certain flux range and 
the cosmic infrared background level. Fig. [I| to Fig. || are 
the model predictions for the European Large Area ISO 
Survey (ELAIS). This survey covered 12 deg 2 of the sky in 
four main areas and was carried out with the ISOPHOT 
instrument onboard the Infrared Space Observatory ISO, 
which is at least an order of magnitude deeper than the 
IRAS 100 pm survey. It therefore provides an important 
constraint for our model of galaxy evolution. The major¬ 
ity of the optical identification of the detected sources are 
for interacting pairs or small groups of galaxies, which 
may indicate that the ELAIS sample includes a significant 
fraction of luminous infrared galaxies from galaxy merg¬ 
ers. Although there is some discrepancy in the data re¬ 


duction, previous estimations show that the source counts 
are mostly in agreement with strongly evolving starburst 
models, with a rapid increase in the fraction of ULIGs 
towards high redshift (Efstathiou et al. 2000, Matsuhara 
et al. 2000, Serjeant et al. 2001). From our calculations, 
we see in Fig. to Fig. that the differential number 
counts of 90 pm, 15 pm and 170 pm for a reliable subset 
of the detected sources could be sufficiently accounted for 
by the infrared burst phase when a population of ultralu¬ 
minous infrared sources with L ir > 10 12 Lq could be pro¬ 
duced by the merger-triggered starburst/AGN activities 
at z ~ 1 (Kawara et al. 1998, Elbaz et al. 1999, Efstathiou 
et al. 2000, Dole et al. 2001). The enormous infrared emis¬ 
sion, especially at far-infrared wavelengths is modelled by 
a modified black body spectrum which we discussed in 
the previous section. We assume that the starburst merg¬ 
ing system has a similar MIR emission feature as Arp220, 
but modified by the observed flux correlation of Sis/Seo 
from IRAS and ISOCAM deep surveys. We found from the 
calculation that the dust temperature of these starburst 
merging systems would be higher than that of the nearby 
starburst ULIGs Arp220, with dust temperature T = 65 K 
and (3 = 1.5 for a best fit result. Fig. [| shows our model 
fitting for the differential counts of the IRAS 60 pm deep 
survey, and Fig. || shows the integrated number counts of 
submillimeter SCUBA deep survey at 850 pm (Hacking et 
al. 1987, Moshir et al. 1992, Barger et al. 1999, Blain et 
al. 1999). Almost all the number counts could be repro¬ 
duced quite well by such an evolutionary scenario, except 
for the ISOCAM 15 pm differential number counts, where 
our model prediction shows a slight excess at the bright 
part of S 15 ~ 2 mJy. The reason could be that we simply 
adopt the mid-infrared emission feature of Arp 220 for the 
case of the starburst merging system around z ~ 1. We 
hope to improve the current results by further theoreti¬ 
cal modelling and additional observational constraints for 
the emission properties at mid infrared bands from future 
infrared missions. 

The infrared background in this calculation gives 
2AnW m~ 2 sr -1 at 15 pm, 1.9 nW m~ 2 sr^ 1 at 60pm, 
3.8nlU?n -2 sr _1 at 90 pm, 10.6 nWm~ 2 sr^ 1 at 170 pm, 
which are all consistent with current upper limits from 
TeV detections, COBE results and the resolved fraction 
of the CIRB by the deep ISO surveys (Funk et al., 1998, 
Guy et al. 2000, Hauser & Dwek 2001). 

The redshift distribution of the ISOCAM 15 pm 
contributing sources within the detected flux range 
(0.1 mJy ~ 10 mJy) from our model calculation is shown 
in Fig. |(]. It gives a rough statistical finding that these 
luminous infrared sources cover a wide redshift range of 
0.5 ~ 2.5, peaking at z ~ 1. Comparing our model pre¬ 
diction and the redshift distribution of 15 pm sources with 
S 15 > 120 pJy in the HDF North and the z-distribution 
of sources in the CFRS field (Flores et al. 1999, Cohen et 
al. 2000, Aussel et al. 2001), we found that the starburst 
mergers at z ~ 1 in our model are good candidates for 
a strongly evolving population that results in the strong 
evolution in mid- and far-infrared deep surveys. Recent 
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redshift estimation from sub-mm follow up of 10 known 
FIRBACK 170gm ISO sources by Scott et al. (2000) sug¬ 
gests that they are in a redshift range of 0 ~ 1.5, still 
consistent with our current model predictions. However, 
these redshift determinations strongly depend on the as¬ 
sumption of the dust properties. We need further accurate 
measurements for the robust constraints of the models. 
We discussed in a previous paper that shifting the peak 
redshift of these ULIGs by a factor of 2 could affect the 
source count fitting of the IRAS 60deep survey, espe¬ 
cially for a low redshift peak (z < 0.5). Strong evolution of 
the ULIGs to z ~ 1 may be the most reasonable case for 
the existing model constraints from both the infrared deep 
surveys and the cosmic infrared background upper limits 
from high energy TeV detections, as well as the indicated 
star formation history by UV/optical deep surveys (Lilly 
et al. 1996, Connolly et al. 1997, Madau et al. 1998). 

We plot the redshift distribution of ULIGs ( vL v > 
10 12 L© ) to understand the evolutionary properties of the 
ULIGs from mergers in our model. A rapid increase in the 
number density of ULIGs up to z ~ 1 is seen in Fig. [t|, 
which is actually consistent with a scenario where galaxy 
merger rates increase dramatically during that epoch as 
seen from various observations and theoretical consider¬ 
ations (Zepf & Koo 1989, Carlberg 1992, Burkey et al. 
1994, Carlberg et al. 1994). However, the number density 
of ULIGs decreases beyond z ~ 2 — 3, which may reflect 
a stage when merger pairs are mostly dwarfs and the in¬ 
frared emissions are less than 10 12 Lq even with intensive 
starburst activities triggered by mergers. In this scenario, 
an infrared luminous tail of the luminosity function may 
form at z ~ 1 , with enormous infrared emission enhance¬ 
ment. 

There is still no firm statistical basis for the classifica¬ 
tion of starbursts and AGNs from current spectroscopies. 
We thus adopt the observed AGN local luminosity func¬ 
tion of Rush et al. (1993) as a model constraint, and as¬ 
sume in our calculations that the observed starburst galax¬ 
ies and Seyferts follow the same evolutionary track, based 
on naive thinking that starburst/AGN may both be trig¬ 
gered by galaxy interactions. We know the subtle differ¬ 
ences in the dust emission properties could result in a 
different fraction of their contribution. This is still far too 
early to discuss here. We thus adopt only the SED of the 
Cloverleaf quasar which represents a phase poor in cold 
gas, as well as the dust enshrouded phase of F10214+4724 
as two typical AGN templates in our calculation. We know 
from the result that the AGN contribution is only a small 
fraction of the whole and our current model prediction 
is within the present understanding of this issue, i.e. the 
starburst powered ULIGs are dominating over the AGN 
powered ones (Fig. [?]) and may take over at higher red- 
shifts and in the higher luminosity case (Lutz et al. 1998, 
Tran et al. 2001). 



S_{ 15} (mJy) 

Fig. 1 . The model prediction of the differential number 
counts of ISOCAM 15 fim normalized to the Euclidean 
law ( dN/dS * S ' 2 5 ). The data points are the normalized 
counts from a variety of ISO deep surveys (Elbaz et al. 
1999). The line represents the sum of the contribution 
from three populations (starburst galaxies, spiral galaxies 
and Seyferts). In our model, we assume that the popula¬ 
tion of starburst galaxies, especially with spheroidal mor¬ 
phology, are the product of galaxy mergers which would 
experience infrared emission enhancement because of the 
merger-triggered starburst activities. Their contribution 
to the ISOCAM 15 /ira deep survey is shown by the dot- 
dashed line from our Monte-Carlo simulation. The dotted 
line corresponds to the non-evolving spiral galaxies, and 
the short dashed line is from Seyferts, which are assumed 
to have the same evolutionary track as starburst galaxies 
in a galaxy evolutionary scheme with galaxy mergers. 

4. Summary 

We have described a galaxy evolutionary scenario with 
galaxy mergers in a CDM cosmology, where starburst and 
AGN activities may be triggered by the merger events. 
With a reasonable assumption of the ultraluminous in¬ 
frared burst phase from gas rich mergers at high redshift, 
we sucessfully interpreted the strong evolution of IRAS 
60 /xm deep survey, leaving the infrared background still 
with a low limit of 1.9 nW m~ 2 sr” 1 , consistent with the 
upper limits from recent TeV 7 ray detection of nearby 
Blazars. 

Lacking a comprehensive understanding of the gas and 
dust environment of faint ISO sources, especially the star- 
burst merging system at z ~ 1 , we adopt the template 
spectral energy distribution such as that of Arp220 by 
Silva et al. (1998) as typical for nearby starburst galax¬ 
ies. Meanwhile, we construct a simple SED for the star- 
burst mergers at z ~ 1. The far-infrared emission in such 
a system is modelled by a single temperature, optically 
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Fig. 2. The fitting of the ELAIS differential source count 
at 90 gm. The data are from C-90 filter of the C100 
ISOPHOT detector array [filled squares: Efstathiou et al. 
2000; open trianguler: Linden-V(/)rnle et al. (2000); star: 
Juvela et al. (2000)]. The meaning of the lines is the same 
as in Fig-S 



Fig. 3. The result of FIRBACK 170 fim ISO deep survey 
differential number count fitting from our model calcula¬ 
tion. The data are from Dole et al. (2001). 

thin dust law with modified black body emission. The 
MIR emission feature is assumed to be similar to Arp 
220, but modified by the flux correlation from IRAS and 
ISOCAM observations. In this case, we can further inves¬ 
tigate such a merger-driven galaxy evolutionary scenario 
at other infrared and submillimeter wavelengths with ISO 
and SCUBA deep surveys. Our calculation shows that 



Fig. 4. The fitting of the IRAS 60 /im source counts from 
three major infrared emitters (starburst galxies, spiral 
galaxies and Seyferts). The source counts of starburst 
galaxies and Seyferts are from the Monte-Carlo simula¬ 
tion where the evolution of both activities are triggered by 
galaxy-galaxy interactions/mergers during structure for¬ 
mation. The spiral galaxy is assumed to have a mild con¬ 
stant star formation history, i.e. a non-evolving popula¬ 
tion in our calculation. The data are from the IRAS Point 
Source Catalogue (1985) (PSC), Hacking et al. IRAS deep 
survey (HCH), FSC from deep surveys by Moshir et al. 
(1992) and Saunders (1990). 



Fig. 5. Integral number counts at 850 /im. The open 
squares are from Blain et al. (1999), and the filled squares 
are from Barger, Cowie and Sanders (1999). The meaning 
of lines are the same as in previous figures. 
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0.1mJy<S_{15}<10mJy 



Fig. 6. The redshift distribution of three infrared contrib¬ 
utors (starburst galaxies, spiral galaxies and AGNs) at a 
flux range of 0.1 mJy ~ 10 mJy at 15 [im from our model 
calculation. The redshift of these ISO far-infrared sources 
would cover a wide range and peak near z ~ 1. We in¬ 
clude also the AGN contribution in our model based on 
the observed AGN Local Luminosity Function from Rush 
et al. (1993). 


the current results of multi-wavelength deep surveys at 
ISOCAM 15 fim, ELAIS 90^m, FIRBACK 170/im, IRAS 
60gm and SCUBA 850 ym number counts could be suf¬ 
ficiently accounted for by the merger-triggered infrared 
enhancement at £ ~ 1 from our model with the dust tem¬ 
perature (T ~ 65 K , and f3 ~ 1.5), slightly higher than the 
local starburst galaxy Arp220. Future accurate redshift 
measurements and multiband photometries would provide 
a robust model check. 

The background levels at these wavelengths are 
estimated from our model calculation, which gives 
2.4nlkm" 2 sr _1 at 15gm, 3.8 nW m~ 2 sr^ 1 at 90 yrri, 
10.6nIU to ” 2 sr ” 1 at 170 /jto, still compatible with the cos¬ 
mic infrared background level both from the upper limit 
of high energy TeV 7 ray detection of nearby Blazars and 
from COBE and ISO results. 

The redshift distribution of the luminous infrared 
sources within the ISOCAM 15 y,m detection flux range 
(0.1 mJy ~ 10 mJy) from our calculation is plotted in 
Fig. 0 The redshift distribution of these sources cover a 
wide redshift range from 0.5 ~ 2.5 and peak around a 
mean redshift of z r-~j \ _ 

We also plot the redshift distribution of ULIGs ( vL v > 
10 12 Lq) in Fig. 0 It shows a strong increase in the ultra- 
luminous infrared population untill a mean redshift z ~ 1 , 
and decreases by a factor of about 2 by 2 ~ 2 — 3. This is 
probably the major difference between our current calcu¬ 
lation and other models, and indicates that the infrared 


Lir> 10 A 12 Lsun 



z 


Fig. 7. The redshift distribution of ultraluminous infrared 
sources (starbursts or AGNs) with Li r > 10 12 Lq from our 
model calculation. The number density of the ultralumi¬ 
nous infrared sources increase dramatically untill z ~ 1 , 
but quickly decrease afterward. 


luminous tail may be produced at the cosmic epoch of 
z ~ 1 , when the merger rate and the size of parent galax¬ 
ies are suitable for such an infrared emission enhancement. 

A brief discussion about the fraction of contribution 
from AGNs and starbursts from our calculation is given 
in section |j. We assumed in the model a similar evolu¬ 
tionary track for the starburst galaxies and AGNs based 
on the idea that both AGNs and starbursts may be trig¬ 
gered by galaxy interactions, where the AGN population 
is constrained by the observed Local Luminosity Function 
of Seyferts from Rush et al. (1993). Given the uncer¬ 
tainty of the dust properties of ULIGs, especially for 
those harboring an AGN in the center, we adopt here 
only two typical SED templates of Cloverleaf QSO and 
F10214+4724. In this case, we give a rough estimation 
of the relative abundance of AGN and starburst-powered 
ULIGs ( L ir > 10 12 Lq) of ~ 1/5, which seems to be close 
to the recent submillimeter observations of Chandra X-ray 
sources (Almaini et al. 1999, Barger et al. 2001, Gunn & 
Shanks 2001). 
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